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Abstract 

We present the microscopic equation of growing interface with quenched noise for 
the Tang and Leschhorn model [Phys. Rev. A 45, R8309 (1992)]. The evolution 
equations for the mean heigth and the roughness are reached in a simple way. 
Also, an equation for the interface activity density (i.e. interface density of free 
sites) as function of time is obtained. The microscopic equation allows us to express 
these equations in two contributions: the diffusion and the substratum one. All the 
equation shows the strong interplay between both contributions in the dynamics. 
A macroscopic evolution equation for the roughness is presented for this model for 
the critical pressure p = 0.461. The dynamical exponent (3 = 0.629 is analytically 
obtained in a simple way. Theoretical results are in excellent agreement with the 
Monte Carlo simulation. 



The growth of interfaces under nonequilibrium conditions is an interesting 
natural phenomenon. The interfaces has been characterized through scaling 
of the interface width w with time t and lateral size L. It is known that w ~ L a 
for t L 01 !^ and w ~ t@ for t <C L a ^. (3 and a are called the dynamical and 
the roughness exponents, respectively. 

The knowledge of the effects of the nonlinearities and the disorder of the media 
are important to describe the growing interface motion. Some experiments 
such as the growth of bacterial colonies and the motion of liquids in porous 
media, where the disorder is quenched, i. e. disorder due to the inhomogeneity 
of the media where the moving phase is propagating, are well described by 
directed percolation depinning (DPD) model proposed simultaneously by Tang 
and Leschhorn (TL) [1] and Buldyrev et. al [2]. In this work we will focus on 
TL model. This model predicts a dynamical exponent (3 ~ 0.63 for q c = 
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1 — p c = 0.539, where p c is the critical pressure. However, some authors [2-4] 
found a dynamical exponent (3 faraway from the criticality. Yang and Hu [4] 
found (5 as function of q, analyzing the roughness as function of t. In a previous 
work, Braunstein and Buceta [5] show that the scaling law holds only at the 
criticality. 

We present now the Microscopic Equation (ME) for the TL model. The inter- 
face growth takes place in a lattice of edge L with periodic boundary condi- 
tions. A random pinning force g(r) uniformly distributed in [0, 1] is assigned 
to every cell of the lattice. For a given pressure p, the cells are divided in 
two groups, active (free) cells with g(r) < p and inactive (blocked) cells with 
g(r) > p. The interface is specified by a set of integer column heights hi 
(i — 1, . . . , L). During the growth, a column is selected at random with prob- 
ability 1/L and compared his height with those of his neighbors. The time 
evolution equation for the interface in a time step St = 1/L is 

hi(t + St) = Kit) + St [W i+1 + Wi-i + F^) Wi] , (1) 

with 

Wi = l- Q{hi - min(/ij_!, h i+1 ) - 2) , 
W i±1 = Q(h i±1 - min(/ii, h i±2 ) - 2){[1 - G(hi - h i±2 )\ + \S hi:hi±2 } . 

Here h'^ = hi + 1 and Q(x) is the unit step function defined as 6(x) = 1 for 
x > and equals to otherwise. F^h'i) equals to 1 if the cell at the height h\ 
is free or active (i.e. the growth may occur at the next step) or if the cell 
is blocked or inactive. Fi is the interface activity function. Taking the limit 
St — > and averaging over the lattice we obtain {h = (hi)) 

^ = (1 - Wi) + (FiWi) . (2) 

This equation allow us the identification of two separate contributions: the 
difussion (1 — Wi) and the substratum one (Fi Wi). Yang and Hu [3] defined 
two kind of growth events: the event in which the growth occurs at the chosen 
site as a type A (defined by us as local growth) and the event in which the 
growth occurs at the adjoint site as type B (our growth by difussion). They 
count, in numerical simulation, the events number N^it) of type A and N B (t) 
of type B in a time interval L. They do not identify this terms as the substra- 
tum and the diffusion contributions to the mean height derivative. Remark 
that N A (t) oc (FiWi) and N B (t) oc (1 - W t ). They concluded that events A 
roughen the interface while events B flatten it. We shall see bellow that this 
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last contribution (diffusion) enhances the roughness near the critical point in 
contradiction with the statement of Yang and Hu. 

We found an amazing numerical result: 

{F i (l-W i ))=p{l-W i ) . (3) 

We could not obtain analytically this result. Using this result, the interface 
activity density (IAD) / = (Fj) is 

f = p{l-W i ) + {F i W i ) . (4) 

The substratum contribution dominates the behavior of / in the earliest 
regime. Figure 1 shows / as function of time. At the initial time, / is equal 
to p and in the earliest regime / decreases as t increases. In the intermediate 
regime (t <C L a ^) only at the critical value a temporal power law holds. 

From Eq. (1), the derivative of the square interface width dw 2 / dt (DSIW) 
can also be expressed by means of two additive contributions. The diffusion 
contribution is 

2 [ ((1 - Wi) minfa-i, h i+1 )) - (1 - Wi) (h,) } , (5) 

and the substratum contribution is 

2[(h i F i W i )-(h i )(F i W i )}. (6) 

At short times, the diffusion process is unimportant because Wi is close to 
one. As t increases, the behavior of this contribution depends on q. Notice, 
from Eq. 5, that the diffusion contribution may be either negative or positive. 
The negative contribution tends to smooth out the surface, this case domi- 
nated for small q. The positive diffusion contribution enhances the roughness. 
This last effect is very important at the critical value. At this value, the sub- 
stratum contribution is practically constant, but the diffusion contribution 
is very strong (see Figure 2), enhancing the roughness in contradiction with 
the statement of Yang and Hu. This last contribution has important duties 
on the power law behavior. Notice that the roughness is the result of the 
competition between both contributions. Generally speaking, the substratum 
roughen the interface while the diffusion flatten it for small q, but the dif- 
fusion also roughen the interface when q increases. The diffusion is enhaced 
by local growth. The growth by diffusion may also increase the probability of 
local growth. This crossing interaction mechanism makes both, local growth 
and growth by diffusion, dominants in turn near the criticality. 
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The ME approach to obtain the Langevin equation for the motion of surface 
has been used in the Das Sarma-Tamborenea model and the ballistic deposi- 
tion with termal noise [6,7]. A perturbative expansion of the ME equation can 
be performed if one regularizes the step function. However, in the TL model 
it is not clear how to treat the quenched noise. The modeling of the quenched 
noise results an open question due to the complex nature of the growth. Never- 
theless, we obtained a phenomenological equation for the roughness temporal 
derivative [5]. The DSIW in the critical regime satisfy 

-^j- = ci exp[-2g(t + *„)] + c 2 (t + t ) 2 ^ 1 , (7) 

where the constants Cj (j = 1,2) and to are related to the dynamics mech- 
anism, to is a small constant (q to <C 1) that takes into account the initial 
condition dw 2 / dt J t =o — P- Introducing the constant c\ = 2 — 1/(3, and using 
the initial condition, at zero order in the exponential expansion around qto = 0, 
we obtain t = [(p — c 1 )/c 2 ] 1 ^ 2/3_1 - ) . The minimum of the DSIW is amazing 
reached at t — tt, as show the Monte Carlo results (see Figure 3). Using the 
second derivative of interface width evaluated in the minimum and replacing 
ci, at zero order in t expansion, we obtain c 2 = (2q/(3) e~ 2g7T ir^ 13 ^. All 
the dynamic process is governed by the competition between the two terms of 
Eq. (7). Figure 3 shows that at t — 1 the contributions are equal. With this 
condition, we obtain the theoretical (3 = 0.629 in excellent agreement with our 
numerical simulation results. The earliest regime is determinated by the two 
terms until t ~ tt. Next the dynamics is governed by the power term. 

Summarizing, we wrote the ME, for the evolution of the height, for the Tang 
and Leschhorn model. The ME allows us to derive the evolution of the mean 
height, the IAD and the DSIW in terms of two separate contributions: the 
substratum and diffusion one. The identification of this two separate contri- 
butions shows the crossing interaction between both. In particular, both are 
important at the criticallity for the roughness. We wrote a phenomenological 
macroscopic equation for the DSIW that allows us to predict the dynamical 
exponent in a simple way. In the present, we are working on the DPD model 
of Buldyrev et. al. Previous results show the same qualitative features for the 
diferents contributions of the ME and for the macroscopic equation as was 
expected for a model that belongs to the same universality class. 
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Fig. 1. j jp versus t. The parameter q is 0.51 (□), 0.539 (O) arm 0.6 (v)- All cases 
shows the same behavior in the early time regime. The subcritical (supercritical) 
case shows that the IAD asymptotically goes to certain constant (zero) . The critical 
case shows that the IAD goes as t~ v (in solid line) with rj = 0.400 ± 10~ 3 . 



Fig. 2. Diffusion contribution to the DSIW (hollow symbols) and the substratum 
one (fill symbols) versus Int for q = 0.539 (circles) and q = 0.3 (triangles). 
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Fig. 3. For the critical q = 0.539 we show the plot of the theorically DSIW versus 
time in solid line, Monte Carlo results in circles and the first (second) term of DSIW 
versus time in dashed line a (b). The vertical solid line is used as a guide to show 
that both contributions are equal at t = 1. The horizontal solid line is used as a 
guide to show the initial condition. 
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